Coulomb blockade of anyons 
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Coulomb interaction turns anyonic quasiparticles of a primary quantum Hall liquid with filling 
factor v = l/(2m + 1) into hard-core anyons. We have developed a model of coherent transport of 
such quasiparticles in systems of multiple antidots by extending the Wigner- Jordan description of 
ID abelian anyons to tunneling problems. We show that the anyonic exchange statistics manifests 
itself in tunneling conductance even in the absence of quasiparticle exchanges. In particular, it can 
be seen as a non- vanishing resonant peak associated with quasiparticle tunneling through a line of 
three antidots. 
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Quasiparticles of two-dimensional (2D) electron liq- 
uids in the regime of the Fractional Quantum Hall effect 
(FQHE) have unusual properties of fractional charge [l[ 
and fractional exchange statistics 0, Q . The fractional 
quasiparticle charge was observed in experiments on an- 
tidot tunneling [j] and shot-noise measurements [B|, @. 
The situation with fractional statistics is so far less cer- 
tain even in the case of the abelian statistics, which is the 
subject of this work. Although the recent experiments Q 
demonstrating unusual flux periodicity of conductance 
of a quasiparticle interferometer can be interpreted as a 
manifestation of the fractional statistics [a, [9|1 , this inter- 
pretation is not universally accepted 

[ml 

111 . There is a 

number of theoretical proposals (see, e.g.. |l2l. [l3j) sug- 
gesting tunnel structures where the statistics manifests 
itself through noise properties. Partly due to complex- 
ity of noise measurements, such experiments have not 
been performed successfully up to now. In this work, 
we show that coherent quasiparticle dynamics in multi- 
antidot structures should provide clear signatures of the 
exchange statistics in dc transport. Most notably, in tun- 
neling through a line of three antidots, fractional statis- 
tics leads to a non- vanishing peak of the tunnel conduc- 
tance which would vanish for integer statistics. 

These effects rely on the ability of quantum anti- 
dots to localize individual quasiparticles of the QH liq- 
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FIG. 1: Tunneling of anyonic quasiparticles between oppo- 
site edges of an FQHE liquid through quasi- ID triple-antidot 
systems: (a) loop, (b) open interval. Quasiparticles tunnel 
between the edges and the antidots with rates Ti^- The an- 
tidots are coupled coherently by tunnel amplitudes A. 
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antidots are very similar to those associated with the 
Coulomb blockade [16| in tunneling of individual elec- 
trons in dots. For instance, similarly to a quantum dot 
[h| , the linear conductance of one antidot shows periodic 
oscillations with each period corresponding to the addi- 
tion of one quasiparticle [3, Q, [H, Q, HH • Recently, we 
have developed a theory of such Coulomb-blockade-type 
tunneling for a double-antidot system [jfij . where quasi- 
particle exchange statistics does not affect the transport. 
The goal of this work is to extend this theory to antidot 
structures where the statistics does affect the conduc- 
tance. The two simplest structures with this property 
consist of three antidots and have quasi-lD geometries 
with either periodic or open boundary conditions (Fig.[T|). 
A technical issue that needed to be resolved to calculate 
the tunnel conductance is that the anyonic field opera- 
tors defined through the Wigner- Jordan transformation 
are not fully sufficient in the situations of 
tunneling. As we show below, to obtain correct matrix 
elements for anyon tunneling, one needs to keep track 
of the appropriate boundary conditions of the wavefunc- 
tions which are not accounted for in the field operators. 

Specifically, we consider the antidots coupled by tun- 
neling among themselves and to two opposite edges of 
the quantum Hall liquid (Fig. [T|) . The edges play the role 
of the quasiparticle reservoirs with the transport voltage 
V applied between them. We assume that the antidot- 
edge coupling is weak and can be treated as a pertur- 
bation. Quasiparticle transport through the antidots is 
governed then by the kinetic equation similar to that for 
Coulomb-blockade transport through quantum dots with 
a discrete energy spectrum [251 ] . Coherent quasiparticle 
dynamics requires that the relaxation rate created by 
direct Coulomb antidot-edge coupling is weak. This con- 
dition should be satisfied if the edge-state confinement 
is sufficiently strong [2(J. The requirement on the con- 
finement is less stringent in the case of the antidot line 
(Fig. lb), in which antidot quasiparticles move along the 
edge, suppressing the antidot-edge coupling at low fre- 
quencies. We also assume that all quasiparticle energies 
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on the antidots, tunnel amplitudes A, temperature T, 
Coulomb interaction energies U between quasiparticles 
on different antidots, are much smaller than the energy 
gap A* for excitations on each antidot. This condition 
ensures that the state of each antidot is characterized 
completely by the occupation number n of its relevant 
quantized energy level. In any given range of the back- 
gate voltage or magnetic field (which produces the over- 
all shift of the antidot energies - see, e.g., [3, EH, EH), 
there can be at most one quasiparticle on each antidot, 
n = 0,1. This "hard-core" property of the quasipar- 
ticles means that they behave as fermions in terms of 
their occupation factors, despite the anyonic exchange 
statistics. All these assumptions can be summarized as: 
T d ,T 3 < A,C/,T< A*. 

Under these conditions, the antidot tunneling is domi- 
nated by the antidot energies. The quasi-lD geometry of 
the antidot systems we consider makes it possible to in- 
troduce the quasiparticle "coordinate" x numbering suc- 
cessive antidots; e.g., x — —1, 0, 1 for systems in Fig. Q] 
The quasiparticle Hamiltonian can be the written as 

H = y^[e x n g - (Ax£ x+ i£x + h.c.)] + ^ U x , y n x n y , (1) 

x x<y 

where e x are the energies of the relevant localized states 
on the antidots (taken relative to the common chemical 
potential of the edges at V = 0), A^ is the tunnel cou- 
pling between them, U XtV is the quasiparticle Coulomb 
repulsion, and n x = Q^ x - The quasiparticle operators 
£| , £e in (fTJ) can be viewed as the Klein factors left in the 
standard operators for the edge-state quasiparticles when 
all the edge magneto-plasmon modes are suppressed by 
the gap A*. Characteristics of such Klein factors depend 
on the geometry of a specific tunne ling problem; non- 
trivial examples can be found in 

[3 El & S3- In the 
Hamiltonian (fT]), £ x describe the hard-core anyons with 
exchange statistics iru. Wigner- Jordan transformation 



expresses them through the Fermi operators c x 2l| : 



c x j £,y£,x — £,x£,y€ 



iTnssgn(x-y) 



(2) 



with similar relations for . 

Anyonic statistics creates an effective interaction be- 
tween the quasiparticles which can be understood as the 
Aharonov-Bohm (AB) interaction between a flux tube 
"attached" to one of the particles and the charge carried 
by another. In general, this interaction can be masked by 
the direct Coulomb interaction U x ,y In the antidot loop 
(Fig. QJl), however, U x . y is constant, U x<y = U, and the 
interaction term in ([T]) reduces to Un(n — l)/2, with n = 
J2 X n x ~ the total number of the quasiparticles on the an- 
tidots. In this case, the Coulomb interaction contributes 
to the energy separation between the group of states with 
different n, but does not affect the level structure for 
given n. The hard-core property of quasiparticles limits n 
to the interval [0, 3] . For n = and n — 3, the system has 



the "empty" and "completely filled" state with respective 
energies E = , E 3 = J2 x e x + 3U. The spectrum En- of 
the three n = 1 states \lk) = Y^ x 4>k{x)£,l\0) , is obtained 
as usual from ([1]). In the uniform case e x = e, A x = A, 
with an external AB phase <p, one has 4>k(x) = e lkx /L 1 ' 2 
and 



Eik = e — A cos k , k — (2irm + ip)/L , 



(3) 



where m = 0, 1, 2, and the loop length is L = 3. 

Anyonic statistics can be seen in the n — 2 states, 
\2l) = (l/V2)E xy Mx,y)4&\°)- The fermion-anyon 
relation ([2]) suggests that the stationary two-quasiparticle 
wavefunctions should coincide up to the exchange phase 
with that for free fermions: 



iTt(l-u)sgn(x-y)/2 
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4>q{x) 4> q (y) 
<f>p(%) <f» P (y) 



(4) 



Here <^s are the single-particle eigenstates of the Hamilto- 
nian (fT]) . (The states ((4]) are numbered with the index I of 
the third "unoccupied" eigenstate of (fT]) complementary 
to the two occupied ones q,p.) The boundary conditions 
for the 0s are affected by the exchange phase in Eq. Q. 
To find them, we temporarily assume for clarity that co- 
ordinates x, y are continuous and lie in the interval [0, L]. 
Subsequent discretization does not change anything sub- 
stantive in this discussion. The ID hard-core particles 
are impenetrable and can be exchanged only by moving 
one of them, say x, around the loop from x = y + to 
x = y — (Fig. [2^i). Since the loop is imbedded in the 
underlying 2D system, such an exchange means that the 
wavefunction acquires the phase factor e I7riy , in which the 
sign of v is fixed by the properties of the 2D system, e.g. 
the direction of magnetic field in the case of FQHE liq- 
uid. Next, if the second particle is moved similarly, from 
y — x + 0toy = x — 0, the wavefunction changes in the 
same way, for a total factor e l27ri/ . Equation (0| shows 
that only one of these changes can agree with the ID form 
of the exchange phase. As a result, the wavefunction (|3|) 
satisfies different boundary conditions in x and y: 

ML,y) = M0,yV v , 1>i{x,l) =^(*,o)e^+ 2 ™). (5) 

Conditions (O on the wavefunction (H|) mean that the 
single-particle functions <j) in (U) satisfy the boundary 
condition that correspond to the effective AB phase 
ip' = if + 7r — ny, i.e. the addition of an extra quasiparti- 
cle to the loop changed the AB phase by ir — nu, where 
—ttv comes from the exchange statistics and it from the 
hard-core condition. This gives the energies of the two- 
quasiparticle states (U]) as U + E\ q + Ei p) where, if the 
loop is uniform, the single-particle energies are given by 
Eq. with ip — > ip' . In this case, J^k^ik — 0, and 
the energies En of the two-quasiparticle states can be 
written as: 



Eo 



2e + U- AcosZ, I = (2irm' + tp - irv)/3 , (6) 
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FIG. 2: Exchanges of hard-core anyons on a ID loop: (a) real 
exchanges by transfer along the loop embedded in a 2D sys- 
tem; (b) formal exchanges describing the assumed boundary 
conditions (JSJ) of the wavefunction. 



where m! — 0, 1, 2. 

One of the consequences of this discussion is that the 
sign of v in the ID exchange phases of Eqs. ^ and ((H) 
can be chosen arbitrarily for a given fixed sign of the 
2D exchange phase. Reversing this sign only exchanges 
the character of the boundary conditions (JSJ) between x 
and y. This fact has simple interpretation. Although 
the ID hard-core anyons can not be exchanged directly, 
formally, coordinates x and y in Eq. ^ are independent 
and one needs to define how they move past each other 
at the point x = y. Depending on whether the a;-particle 
moves around y from below or (as in Fig.[2jD) from above, 
its trajectory does or does not encircle the flux carried 
by the y particle, and the boundary condition for x is 
or is not affected by the statistical phase. The choice 
made for x immediately implies the opposite choice for 
y (Fig. [2)), accounting for different boundary conditions 
([5]). This interpretation shows that in calculation of any 
matrix elements, the participating wavefunctions should 
be taken to have the same boundary conditions. While 
this requirement is natural for processes with the same 
number of anyons, it is less evident for tunneling that 
changes the number of anyons in the system. Indeed, the 
most basic, tunnel-Hamiltonian, description of tunneling 
into the point z of the system leads to the states 



el|ifc> = (iA/2)^fc(z,y)4^l°> 



(7) 



^ k {x,y) = [<t> k {x)8 y , z - ^ x - v > 9n{x -^8 x , z (j> k {y)}/V2. 

One can see that Eq. ([7]) automatically implies specific 
choice of the boundary conditions which corresponds to 
the tunneling anyon not being encircled by anyons al- 
ready in the system. This means that in the calculation 
of the tunnel matrix elements with the states ([4]), one 
should always pair the coordinate of the tunneling anyon 
with the discontinuous one in ([5]). Then, the tunnel ma- 
trix elements arc obtained as 



{2l\^l\lk) = ^Y j ri{x,z)cj> k {x). 



(8) 



For instance, in the case of uniform loop with states ([3|) 
and |6]), we get up to an irrelevant phase factor 



(2J|£t|ifc) = (2/3)cos[(fc-0/2]. 



(9) 



Specific anyonic interaction between quasiparticles can 
be seen in the fact that the matrix elements © do not 



vanish for any pair of indices fc, I. In the fermionic case 
v = 1, one of the elements © always vanishes for any 
given fc, since the two-particle state after tunneling neces- 
sarily has one particle in the original single-particle state. 
By contrast, the tunneling anyon can shift existing par- 
ticle out of its state. 

The matrix elements involving empty or fully occu- 
pied states coincide with those for fermions. Taken to- 
gether with Eqs. ((8|) and (J9j) for transitions between the 
partially filled states, they determine the rates Tj(E) = 
^jfv (E)\ (£]) | 2 of tunneling between the jth edge and the 
antidots, where is the overall magnitude of the tun- 
neling rate, and 

f v (E) = {2TrT/uj c y- 1 \T{v/2 + iE/2iiT)\ 2 e- E / 2T /2nT{v) 

is its energy dependence associated with the Luttinger- 
liquid correlations in the edges [28(. Here T(z) is the 
gamma-function and io c is the cut-off energy of the edge 
excitations. The rates IT 'j(E) can be used in the stan- 
dard kinetic equation to calculate the conductance of the 
antidot system [2(|. Anyonic statistics of quasiparticles 
affects the position and amplitude of the conductance 
peaks through the shift of the energy levels by quasiparti- 
cle tunneling (described, e.g., by Eq. (J6])) and through the 
kinetic effects caused by the anyonic features in the ma- 
trix elements ([8]). In the case of the antidot loop (Fig.QJ,), 
however, effects of statistics are masked by the external 
AB flux (p through the loop. Since the area of practical 
antidots is much larger than the internal area of the loop, 
ip is essentially random and can not be controlled by ex- 
ternal magnetic field on the relevant scale of one period 
of conductance oscillations. Below, we present the results 
for conductance for the similar case of a line of antidots 
(Fig. [TJa), the conductance of which is insensitive to the 
AB flux, and shows effects of fractional statistics in the 
tunneling matrix elements. 

As before, the quasiparticle Hamiltonian is given by 
Eq. Ilj. In this geometry, the interaction energy U\ = 
Ui t o = J7o,-i between the nearest-neighbor antidots is in 
general different from the interaction Ui = Ux,-i between 
the quasiparticles at the ends. The localization energies 
on the antidots can be written as ej = e + xS + 2X\x\. We 
consider first the unbiased line, (5 = 0. At low tempera- 
tures, T <C A, U, only the ground states of n quasipar- 
ticles with energies E n participate in transport: Eq = 0, 
Ex = e + A - u, E 2 = 2e + 3A - Q + (U a + U b )/2, and 
E 3 = 3e + 2U a + U b + 4A, where w = (A 2 + Aj + A 2 ) 1 / 2 
and uj is given by the same expression with A replaced by 
A = A — (Ux — U%)l2, In this regime, the linear conduc- 
tance G consists of three peaks, with each peak associated 
with addition of one more quasiparticle to the antidots, 



G 



(ev) 2 7172 



o n fv{E n +x — E n ) 



(10) 



T 7i + 72 1 + exp[-(E n+1 - E n )/T] 

where a n = \(n + l|£g|n)| 2 . The amplitudes do, a 2 are 
effectively single-particle, and thus, independent of the 
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FIG. 3: Linear conductance G of the antidot line in a v = 1/3 
FQHE liquid (Fig. [lb) as a function of the common anti- 
dot energy e relative to the edges. In contrast to electrons 
(y = 1, left inset), tunneling of quasiparticles with fractional 
exchange statistics produces non- vanishing conductance peak 
associated with transition between the ground states of one 
and two quasiparticles. The maximum of this peak is shown in 
the right inset (y = 1/3 - solid, v = 1 - dashed line) as a func- 
tion of the bias <5. The curves are plotted for Ai = A2, A = 0, 
7i = 72; conductance is normalized to Go = (e^) 2 ri(0)/Ai. 



exchange statistics: ao = (oj + A)/2w, and 02 = (Co — 
\)/2uj. By contrast, the amplitude a\ of the transition 
from one to two quasiparticles is multi-particle, and is 
found from Eqs. (j4]) and (jHJ) to be strongly statistics- 
dependent, 



a 1 



A 2 A 2 
12 



(uj + \)uj(uj — A)ct> 



cos 2 (ttv/2) 



(11) 



In particular, a\ vanishes for electron tunneling (y = 1), 
but is non-vanishing in the case of fractional statistics, 
e.g., for v = 1/3, when cos 2 (7w/2) = 3/4. This is illus- 
trated in Fig. [3] which shows the conductance G obtained 
by direct solution of the full kinetic equation for tunnel- 
ing through the antidots. Qualitatively, the vanishing 
amplitude a\ for electrons can be understood as a re- 
sult of destructive interference between the two terms in 
the two-particle wavefunction which correspond to dif- 
ferent ordering of the added/existing electron on the an- 
tidot line. The opposite signs of these two terms lead 
to vanishing overlap with the single-particle state in the 
tunnel matrix element. Fractional statistics of quasipar- 
ticles makes this destructive interference incomplete. Fi- 
nite bias <5 7^ along the line suppresses this interference 
making the effect of the statistics smaller. One can still 
distinguish the fractional statistics by looking at the de- 
pendence of the amplitude of the middle peak of conduc- 
tance on the bias 5 shown in the right inset in Fig. [3l 

In conclusion, we have developed a model of coherent 
transport of anyonic quasiparticles in systems of mul- 
tiple antidots. In antidot loops, addition of individual 
quasiparticles shifts the quasiparticle energy spectrum by 
adding statistical flux to the loop. In the case without 
loops, energy levels are insensitive to quasiparticle statis- 
tics, but the statistics still manifests itself in the quasi- 
particle tunneling rates and hence dc tunnel conductance 



of the antidot system. 
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